Magnetic Response in Mesoscopic Hubbard Rings: A Mean Field Study 
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The present work proposes an idea to remove the long standing controversy between the calculated 
and measured current amplitudes carried by a small conducting ring upon the application of an 
Aharonov-Bohm (AB) flux (f>. Within a mean field Hartree-Fock (HF) approximation we numerically 
calculate persistent current, Drude weight, low-field magnetic susceptibility and related issues. Our 
analysis may be inspiring for studying magnetic response in nano-scale loop geometries. 
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I. INTRODUCTION 



In mesoscopic range phase coherence of electronic 
states is of fundamental importance and the existence of 
dissipationlcss current in a mesoscopic conducting ring 
threaded by an Aharonov-Bohm (AB) flux (/> is a spec- 
tacular consequence of quantum phase coherence. The 
existence of persistent current in mesoscopic rings has 
been addressed several years ago in the pioneering work 
of Biittiker, Imry and Landauer [l|. Later, many excel- 
lent experiments (2|-|5| have been done in different sys- 
tems to reveal the phenomenon of persistent current. 
Though in literature many theoretical ^ i22,] as well as 
experimental papers [245| on persistent current are avail- 
able, yet lot of controversies are still present between 
the theory and experiment, and the complete knowledge 
about it in this scale is not very well established even 
today. The unexpectedly large amplitudes of measured 
currents provide to the most important controversy. It 
has been proposed that the electron-electron (c-c) inter- 
actions contribute significantly to the average currents. 
An explanation based on the perturbative calculation in 
presence of interaction and disorder has been proposed 
and it seems to give a quantitative estimate closer to the 
experimental results, but still it is less than the measured 
currents by an order of magnitude, and the interaction 
parameter used in the theory is not well understood phys- 
ically. Though an attempt has been made to explain the 
enhancement of current amplitude by some theoretical 
arguments but the sign of low-field currents cannot be 
predicted precisely and it is an important controversial 
issue between theoretical and experimental results. 

Motivated with these open challenges in the present 
paper we address magnetic response in mesoscopic Hub- 
bard rings threaded by AB flux </>. We try to propose 
an idea to remove the unexpected discrepancy between 
the calculated and measured current amplitudes by in- 
corporating the effect of second-neighbor hopping (SNH) 
in addition to the traditional nearest-neighbor hopping 
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(NNH) integral in the tight-binding Hamiltonian. Using 
a generalized Hartree-Fock (HF) approximation |23l - [25j . 
we numerically compute persistent current (/), Drude 
weight (-D) and low-field magnetic susceptibility (x) as 
functions of AB flux (f), total number of electrons N^ and 
system size N . With this (HF) approach one can study 
magnetic response in a much larger system since here 
a many-body Hamiltonian is decoupled into two effec- 
tive one-body Hamiltonians. One is associated with up 
spin electrons and other is related to down spin electrons. 
But the point is that, the results calculated using gener- 
alized HF mean-fleld theory may deviate from exact re- 
sults with the reduction of dimensionality. So we should 

<3> 

FIG. 1: (Color online). Schematic view of a ID mesoscopic 
ring penetrated by a magnetic flux </>. The filled black circles 
correspond to the positions of the atomic sites. 



take care about the mean-field calculation, specially, in 
one-dimensional systems. To trust our predictions, in the 
present work also we make a comparative study between 
the results obtained from mean-field theory and exactly 
diagonalizing the full many-body Hamiltonian. The later 
approach where a complete many-body Hamiltonian is 
diagonalized to get energy eigenvalues is not suitable to 
study magnetic response in larger systems since the size 
of the matrices increases very sharply with the total num- 
ber of up and down spin electrons. Our results can be 
utilized to explore magnetic response in any interacting 
mesoscopic system. 

We organize the paper as follows. Following a brief 
introduction (Section I), in Section II we describe the 
geometric model and generalized Hartree-Fock theory to 
study magnetic response in the model quantum system. 
Section HI contains the numerical results, and finally, 
summary of our work will be available in Section IV. 



II. MODEL AND THEORETICAL 
FORMULATION 

We start by referring to Fig. [1] where a normal metal 
ring is threaded by a magnetic flux (f>. To describe the 
system we use a tight-binding framework. For a A^-site 
ring, penetrated by a magnetic flux (j) (measured in unit 
of the elementary flux quantum 0o = ch/e), the tight- 
binding Hamiltonian in Wannier basis looks in the form, 
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where, e,;(j is the on-site energy of an electron at the site i 
of spin a (t, ^) . The variable t corresponds to the nearest- 
neighbor (j = i ± 1) hopping strength, while ii gives the 
second- neighbor (fc — i±2) hopping integral. 6 = 2tt(I)/N 
and 9i = Ancp/N are the phase factors associated with the 
hopping of an electron from one site to its neighboring 
site and next-neighboring site, respectively, cj^ and Cia 
are the creation and annihilation operators, respectively, 
of an electron at the site i with spin a. U is the strength 
of on-site Hubbard interaction. 



A. Decoupling of the interacting Hamiltonian 

In order to determine the energy eigenvalues of the in- 
teracting model quantum system described by the tight- 
binding Hamiltonian given in Eq. [I] first we decouple 
the interacting Hamiltonian using generalized Hartree- 
Fock approach, the so-called mean field approximation. 
In this approach, the full Hamiltonian is completely de- 
coupled into two parts. One is associated with the up- 
spin electrons, while the other is related to the down- 
spin electrons with their modified site energies. For 
up and down spin Hamiltonians, the modified site en- 
ergies are expressed in the form, eL = e^^ + U{ni{) and 
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operator. With these site energies, the full Hamiltonian 
(Eq. [1} can be written in the decoupled form as, 
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where, i?-|- and H^ correspond to the effective tight- 
binding Hamiltonians for the up and down spin electrons, 
respectively. The last term is a constant term which pro- 
vides an energy shift in the total energy. 



B. Self consistent procedure 

With these decoupled Hamiltonians [H^ and Hi) of up 
and down spin electrons, now we start our self consistent 
procedure considering initial guess values of (jT-if) and 
(nil). For these initial set of values of (rii-f) and (nii), 
we numerically diagonalize the up and down spin Hamil- 
tonians. Then we calculate a new set of values of (rii-^) 
and (nil). These steps are repeated until a self consistent 
solution is achieved. 



C. Calculation of ground state energy 

Using the self consistent solution, the ground state en- 
ergy Eq for a particular filling at absolute zero tempera- 
ture (T = OK) can be determined by taking the sum of 
individual states up to Fermi energy (Ep) for both up 
and down spins. Thus, we can write the final form of 
ground state energy as, 

Eo^Y ^«t + E ^»i - E U{n^t){n^i) (3) 



where, the index n runs for the states upto the Fermi 
level. Enf (Eni) is the single particle energy eigenvalue 
for n-th cigenstate obtained by diagonalizing the Hamil- 
tonian H^ (Hi). 



D. Calculation of persistent current 

At absolute zero temperature, total persistent current 
of the system is obtained from the expression, /(</)) = 
— c dEo{(f))/d(f) where, Eo{(f>) is the ground state energy. 



E. Calculation of Drude weight 

The Drude weight for the ring can be calculated 
through the relation. 
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where, N gives total number of atomic sites in the ring. 
Kohn [261 has shown that for an insulating system D 
decays exponentially to zero, while it becomes finite for 
a conducting system. 



F. Determination of low-field magnetic 
susceptibility 

The general expression of magnetic susceptibility x at 
any flux (p is written in the form, 
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N^ f dl{<j>) \ 
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(5) 



Evaluating the sign of x((/)) we can able to predict 
whether the current is paramagnetic or diamagnetic in 
nature. Here we will determine x(0) only in the limit 
(/)—>■ 0, since we arc interested to know the magnetic 
response in the low- field limit. 

In the present work we perform all the essential fea- 
tures of magnetic response at absolute zero temperature 
and use the units where c = h ~ e ~ 1. Throughout 
our numerical work we set the nearest-neighbor hopping 
strength t = —1 and second-neighbor hopping strength 
ii = —0.7. Energy scale is measured in unit of t. 



III. NUMERICAL RESULTS AND DISCUSSION 

A. Perfect Hubbard Rings Described with NNH 
Integral 

For perfect rings we choose £,-[- = fij. = for all i and 
since here we consider the rings with only NNH integral, 
the second-neighbor hopping strength ti is fixed at zero. 



1. Energy-flux characteristics 

To explain the relevant features of magnetic response 
we begin with the energy-flux characteristics. As illustra- 
tive examples, in Fig. [2] we plot the ground state energy 
levels as a function of magnetic flux cj) for some typical 
mesoscopic rings in the half-filled case, where (a) and 
(b) correspond to A^ = 5 and 6, respectively. The red 
curves represent the energy levels for the non-interacting 
{U ~ 0) rings, while the green and blue lines correspond 
to the energy levels for the interacting rings where the 
electronic correlation strength U is fixed to 1 and 2, 
respectively. From the spectra it is observed that the 
ground state energy level shifts towards the positive en- 
ergy and it becomes much fiatter with the increase of 
the correlation strength U . Both for the two different 
ring sizes {N = 5 and 6) the ground state energy levels 
vary periodically with AB fiux 0, but a significant differ- 
ence is observed in their periodicities depending on the 
oddness and evenness of the ring size A^. For A^ = 6 
(even), the energy levels show conventional (J)q (= 1, in 
our chosen unit system c = e = h = 1) flux-quantum 
periodicity. On the other hand, the period becomes half 
i.e., 00/2 for TV = 5 (odd). This 0o/2 periodicity disap- 
pears as long as the filling is considered away from the 




half-filling. At the same time, it also vanishes if impu- 
rities are introduced in the system, even if the ring is 
half-filled with odd A^. Therefore, 00/2 periodicity is a 
special feature for odd half-filled perfect rings irrespec- 
tive of the Hubbard strength U , while for all other cases 
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FIG. 2: (Color online). Ground state energy levels as a func- 
tion of flux <j) for some typical mesoscopic rings in half-filled 
case. The red, green and blue curves correspond to [7 = 0, 1 
and 2, respectively, (a) A'' = 5 and (b) A'' = 6. 

traditional 0o periodicity is obtained. 

To judge the accuracy of the mean-field calculations 
in our ring geometry, in Fig. [3] we show the variation of 
lowest energy levels where the eigenenergies are deter- 
mined through exact diagonalization of the full many- 
body Haniiltonian for the identical rings as given in 
Fig. [21 considering the same parameter values. Com- 
paring the results presented in Figs. [5] and [31 we see that 
the mean-field results agree very well with the exact di- 
agonalization method. Thus we can safely use mean-field 
approach to study magnetic response in our geometry. 



2. Current-flux characteristics 

Following the above energy-flux characteristics now we 
describe the behavior of persistent current in mesoscopic 
Hubbard rings. As representative examples, in Fig. [31 we 
display the variation of persistent currents as a function 
of flux (j) for some typical single-channel mesoscopic rings 
in the half-filled case, where (a) and (b) correspond to 
A^ = 15 and 20, respectively. The red, green and blue 
curves in Fig. [4l^a) correspond to the currents for t/ = 0, 
1.5 and 2, respectively, while these curves in Fig. [IJb) 
represent the currents for C/ = 0, 1 and 1.5, respectively. 
In the absence of any e-e interaction {U = 0), persistent 



current shows saw-tooth hke nature as a function of flux 
(f) with sharp transitions at n(j)Q/'2. (red hue of Fig. IH^a)) 
or 7100 (red hne of Fig. IDJb)), where n being an integer, 
depending on whether TV is odd or even. The saw-tooth 
Hke behavior disappears as long as the electronic cor- 
relation is introduced into the system. This is clearly 
observed from the green and blue curves of Fig. SI Ad- 
ditionally, in the presence of C/, the current amplitude 
gets suppressed compared to the current amplitude in 
the non-interacting case, and it decreases gradually with 
increasing U . This provides the lowering of electron mo- 
bility with the rise of U and the reason behind this can 
be much better understood from our forthcoming discus- 
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FIG. 3: (Color online). Ground state energy levels as a func- 
tion of flux (f) for some typical mesoscopic rings in half-filled 
case, where eigenenergies are determined through exact diago- 
nalization of the full many-body Hamiltonian. The red, green 
and blue curves correspond to (7 = 0, 1 and 2, respectively. 
(a) iV = 5 and (b) A^ = 6. 

sion. Both for two different rings with sizes A^ = 15 (odd) 
and 20 (even), persistent currents vary periodically with 
AB flux (/) showing different periodicities, following the 
energy-flux characteristics. For A^ ~ 15, current shows 
00 /2 flux-quantum periodicity, while for the other case 
{N = 20), current exhibits 00 flux-quantum periodicity. 



3. Variation of electronic mobility: Drude weight 

To reveal the conducting properties of Hubbard rings, 
we study the variation of Drude weight D for these sys- 
tems. Drude weight can be calculated by using Eq. S) 
Finite value of D predicts the metallic phase, while for 
the insulating phase it drops exponentially to zero [26| . 

As illustrative examples, in Fig. [S] we show the varia- 
tion of Drude weight I? as a function of electronic correla- 




tion strength U for some typical single- channel Hubbard 
rings. In Fig.ISja) the results are shown for three different 
half-filled rings, where the red, green and blue lines corre- 
spond to the rings with N = 10, 30 and 50, respectively. 
From the curves it is evident that for smaller values of C/, 
the half-filled rings show finite value of D which reveals 
that they arc in the metallic phase. On the other hand, 
D drops sharply to zero when U becomes high. Thus 
the rings become insulating when U is quite large. The 
results for the non-half filled case are shown in Fig. [5l^b), 
where we fix the ring size TV = 20 and vary the electron 
filling. The red, green and blue curves represent TVg = 10, 
14 and 18, respectively, where TVg gives the total number 
of electrons in the ring. For these three choices of TVg, 
the ring is always less than half-filled (since TVg < TV) 
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FIG. 4: (Color online). Persistent current as a function of 
fiux (j) for single-channel mesoscopic rings in half-filled case, 
(a) TV = 15. The red, green and blue curves correspond to 
f/ = 0, 1.5 and 2, respectively, (b) TV = 20. The red, green 
and blue curves correspond to [/ = 0, 1 and 1.5, respectively. 

and the ring is in the conducting phase irrespective of 
the correlation strength U. Now we try to justify the 
dependence of the Hubbard strength U on the electronic 
mobility for these different fillings. To understand the 
effect of U on electron mobility here we measure a quan- 
tity called 'average spin density' (ASD) which is defined 
by the factor ^,- Kn^-^ — ni^)\/N . The integer i is the site 
index and it runs from 1 to TV. By calculating ASD we 
can estimate the occupation probability of electrons in 
the ring and it supports us to explain whether the ring 
lies in the metallic phase or in the insulating one. For 
the rings those are below half-filled, ASD is always less 
than unity irrespective of the value of U as shown by 
the curves in Fig.ISJb). It reveals that for these systems, 
ground state always supports an empty site and elec- 
tron can move along the ring avoiding double occupancy 



of two different spin electrons at any site i in the pres- 
ence of e-e correlation which provides the metallic phase 
{D > 0). For a fixed ring size and a particular strength 
of t/, the ASD increases as the filling is increased towards 
half-filling which is noticed by comparing the three dif- 
ferent curves in Fig. Mjo). On the other hand, in the 
half-filled rings, ASD is less than unity for small value of 
[/, while it reaches to unity when U is large. This be- 
havior is clearly shown by the curves given in Fig. |6l[a), 
where the red, green and blue lines correspond to ASDs 
for the half-filled rings with N = 10, 30 and 50, respec- 
tively. Thus, for low U there is some finite probability of 
getting two opposite spin electrons in a same site which 
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FIG. 5: (Color online). Drude weight as a function of Hub- 
bard interaction strength U for single-channel mesoscopic 
rings, (a) Half-filled case. The red, green and blue curves 
correspond to A'^ = 10, 30 and 50, respectively, (b) Non- 
half-filled case with A'^ = 20. The red, green and blue curves 
correspond to A''e — 10, 14 and 18, respectively. 

allows electrons to move along the ring and the metallic 
phase is obtained. But for large U, ASD reaches to unity 
which means that each site is singly occupied either by 
an up or down spin electron with probability 1. In this 
case ground state does not support any empty site and 
due to strong repulsive e-e correlation one electron sit- 
ting in a site does not allow to come other electron with 
opposite spin from the neighboring site which provides 
the insulating phase {D = 0). The situation is somewhat 
analogous to Mott localization in one-dimensional infinite 
lattices. In perfect Hubbard rings the conducting nature 
has been studied exactly quite a long ago using the ansatz 
of Bethe by Shastry and Sutherland [23| . They have cal- 
culated charge stiffness constant (Dc) and have predicted 
that Dc goes to zero as the system approaches towards 
half-filling for any non-zero value of U. Our numerical 
results clearly justify their findings. 



4- Low-field magnetic susceptibility 

Now, we discuss the variation of low-field magnetic sus- 
ceptibility which can be calculated from Eq. [5]by setting 
(^ — 7- 0. With the help of this parameter we can jus- 
tify whether the current is paramagnetic (-|-ve slope) or 
diamagnetic (— ve slope) in nature. For our illustrative 
purposes, in Fig. [7] we show the variation of low-field 
magnetic susceptibility with system size N for some typ- 
ical single-channel mesoscopic rings in the half-filled case. 
Figure [TJ^a) correspond to the variation of low-field mag- 
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FIG. 6: (Color online). Average spin density (ASD) as a 
function of Hubbard interaction strength U for single-channel 
mesoscopic rings, (a) Half-filled case. The red, green and 
blue curves correspond to A'^ = 10, 30 and 50, respectively, 
(b) Non-half-filled case with N = 20. The red, green and 
blue curves correspond to N^ = 10, 14 and 18, respectively. 

netic susceptibility for the non-interacting ([/ ~ 0) rings, 
where the ring size can by anything i.e., either odd, fol- 
lowing the relation N ~ 2n-\-l (n is an integer) , or even, 
obeying the expression N ~ 2n + 2. It is observed that 
both for odd and even N, low- field current exhibits dia- 
magnetic nature. The behavior of the low-field currents 
changes significantly when the e-e interaction is taken 
into account. Depending on the ring size iV, the sign 
becomes -l-ve and — ve as shown by the curves given in 
Figs. [7Jb)-(d). For the interacting rings where the rela- 
tion A^ = 4n-f 2 is satisfied, the low- field current becomes 
diamagnetic (Fig.[7l[b)). The sign becomes paramagnetic 
when N = An (Fig. Eljc)) and A^ == 2?i + 1 (Fig. Uld)). 
Thus, in brief, we say that for non-interacting half-filled 
rings low-field current exhibits diamagnetic response ir- 
respective of A^ i.e., whether A^ is odd or even. For the 
interacting half-filled rings with odd A", low-field current 
provides only the paramagnetic behavior, while for even 
A^, depending on the particular value of A^, the response 



becomes either diamagnetic or paramagnetic. These na- 
tures of low-field currents change for the cases of other 
electron fillings. Hence, it can be emphasized that the be- 
havior of the low-field currents is highly sensitive on the 
Hubbard correlation, electron filling, evenness and odd- 
ness of N, etc. The behavior of zero-field magnetic sus- 
ceptibility in Hubbard rings has been studied extensively 
quite a long back using the Bethe ansatz by Shiba |28| . 
In this work, he has studied magnetic susceptibility per 
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FIG. 7: (Color online). Low-field magnetic susceptibility as a 
function of system size A*' for single-channel mesoscopic rings 
in half-filled case, (a) [7 = 0. A*' is an odd {2n+ 1) or an even 
(2n -I- 2) number, where n is an integer, (b) U — 1. A'' is an 
even number obeying the relation A'' = 4n -|- 2. (c) Z7 = 1. A'' 
is an even number satisfying the relation A'" = 4n. (d) U — 1. 
N is an odd number following the relation A' = 2n + 1. 

electron as functions of electron filling and Hubbard cor- 
relation strength and provided several interesting results. 
From his findings we can clearly justify our presented re- 
sults. 



B. Disordered Hubbard Rings Described with 
NNH and SNH Integrals 

Finally, we explore the combined effect of electron- 
electron correlation and second-neighbor hopping (SNH) 
integral on persistent current in disordered mesoscopic 
rings. 

To get a disordered ring, we choose site energies (e^-i- 
and eii) randomly from a "Box" distribution function of 
width W. As the site energies are chosen randomly it 
is needed to consider the average over a large number of 
disordered configurations (from the stand point of statis- 
tical average). Here, we determine the currents by taking 
the average over 50 random disordered configuration in 



each case to achieve much accurate results. 

As illustrative examples, in Fig.|S]we display the varia- 
tion of persistent currents for some single-channel meso- 
scopic rings considering 1/3 electron filling. In (a) the 
results are given for the rings characterized by the NNH 
integral model. The red curve represents the current for 
the ordered {W — 0) non-interacting (U = 0) ring. It 
shows saw-tooth like nature with AB flux </> providing 0o 
fiux-quantum periodicity. The situation becomes com- 
pletely different when impurities are introduced in the 
ring as clearly seen by the other two colored curves. The 
green curve represents the current for the case only when 
impurities arc considered but the effect of Hubbard inter- 
action is not taken into account. It varies continuously 
with and gets much reduced amplitude, even an or- 
der of magnitude, compared to the perfect case. This is 
due to the localization of the energy eigenstates in the 
presence of impurity, which is the so-called Anderson lo- 
calization. Hence, a large difference exists between the 
current amplitudes of an ordered and disordered non- 
interacting rings and it was the main controversial is- 
sue among the theoretical and experimental predictions. 
Experimental results suggest that the measured current 
amplitude is quite comparable to the theoretically esti- 
mated current amplitude in a perfect system. To remove 
this controversy, as a first attempt, we include the effect 
of Hubbard interaction in the disordered ring described 
by the NNH model. The result is shown by the blue 
curve where U is fixed at 0.5. It is observed that the 
current amplitude gets increased compared to the non- 
interacting disordered ring, though the increment is too 
small. Not only that the enhancement can take place 
only for small values of U, while for large enough U the 
current amplitude rather decreases. This phenomenon 
can be explained as follows. For the non-interacting dis- 
ordered ring the probability of getting two opposite spin 
electrons becomes higher at the atomic sites where the 
site energies are lower than the other sites since the elec- 
trons get pinned at the lower site energies to minimize 
the ground state energy, and this pinning of electrons 
becomes increased with the rise of impurity strength W. 
As a result the mobility of electrons and hence the cur- 
rent amplitude gets reduced with the increase of impurity 
strength W. Now, if we introduce electronic correlation 
in the system then it tries to depin two opposite spin 
electrons those are situated together due to the Coulomb 
repulsion. Therefore, the electronic mobility is enhanced 
which provides quite larger current amplitude. But, for 
large enough interaction strength, mobility of electrons 
gradually decreases due to the strong repulsive interac- 
tion. Accordingly, the current amplitude gradually de- 
creases with U. So, in short, we can say that within the 
nearest-neighbor hopping (NNH) model electron-electron 
interaction does not provide any significant contribution 
to enhance the current amplitude, and hence the contro- 
versy regarding the current amplitude still persists. 

To overcome this controversy, finally we make an at- 
tempt by incorporating the effect of second- neighbor hop- 



ping (SNH) integral in addition to the nearest-neighbor 
hopping (NNH) integral. With this modification a sig- 
nificant change in current amplitude takes place which 
is clearly observed from Fig. [SJb). The red curve refers 
to the current for the perfect {W = 0) non-interacting 
(U ~ 0) ring and it achieves much higher amplitude com- 
pared to the NNH model (see red curve of Fig. [SJa)). 
This additional contribution comes from the SNH inte- 
gral since it allows electrons to hop further. The main 
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of Fig.m^a)). This is due the fact that the SNH integral 
tries to delocalize the electronic states, and therefore, 
the mobility of the electrons is enriched. The situation 
becomes more interesting when we include the effect of 
Hubbard interaction. The behavior of the current in the 
presence of interaction is plotted by the blue curve of 
Fig.[5|b) where we fix C/ = 1.5. Very interestingly we see 
that the current amplitude is enhanced moderately and 
quite comparable to that of the perfect cylinder. There- 
fore, it can be predicted that the presence of SNH integral 
and Hubbard interaction can provide a persistent current 
which may be comparable to the measured current am- 
plitudes. In this presentation we consider the effect of 
only SNH integral in addition to the NNH model, and, 
illustrate how such a higher order hopping integral leads 
an important role on the enhancement of current ampli- 
tude in presence of Hubbard correlation for disordered 
rings. Instead of considering only the SNH integral we 
can also take the contributions from all possible higher 
order hopping integrals with reduced hopping strengths. 
Since the strengths of other higher order hopping inte- 
grals are too small, the contributions from these factors 
arc reasonably small and they will not provide any sig- 
nificant change in the current amplitude. Finally, we can 
say that further studies are needed by incorporating all 
these factors. 



IV. CLOSING REMARKS 



FIG. 8: (Color online). Persistent current as a function of flux 
(jf> for single-channel mesoscopic rings with A'^ = 15 considering 
1/3 electron filling, (a) Rings with only NNH integral. The 
red line corresponds to the ordered non-interacting ring, while 
the green and blue lines correspond to the disordered (W — 2) 
rings with U — and 0.5, respectively, (b) Rings with NNH 
and SNH integrals. The red line represents the ordered non- 
interacting ring, whereas the green and blue line correspond to 
the disordered {W = 2) rings with U = and 1.5, respectively. 

focus of this sub-section is to interpret the combined ef- 
fect of SNH integral and Hubbard correlation on the en- 
hancement of persistent current in disordered ring. To 
do this first we narrate the effect of SNH integral in dis- 
ordered non-interacting ring. The nature of the current 
for this particular case is shown by the green curve of 
Fig. [IJb). It shows that the current amplitude gets re- 
duced compared to the perfect case (red line), which is 
expected, but the reduction of the current amplitude is 
very small than the NNH integral model (see green curve 



To summarize, we have studied magnetic response in 
mesoscopic Hubbard rings threaded by Aharonov-Bohm 
flux <j>. Using the generalized Hartree-Fock (HF) approx- 
imation, we have numerically computed persistent cur- 
rent, Drude weight, low-field magnetic susceptibility and 
some other related issues. Most importantly, we have 
tried to implement an idea to remove the long standing 
problem between the experimentally measured and the- 
oretically predicted current amplitudes. We believe the 
present analysis is found to exhibit several interesting re- 
sults which have so far remained unaddressed. 
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